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Abstract. In a previous work (Mbelek 2001), we modelled the rotation curves (RC) of spiral galaxies by including 
in the equation of motion of the stars the dynamical terms from an external real self-interacting scalar field, if), 
minimally coupled to gravity and which respects the equivalence principle in the weak fields and low velocity 
approximation. This model appeared to have three free parameters : the turnover radius, ro, the maximum tan- 
gential velocity, vg max — vg(ro), plus a strictly positive integer, n. Here, we propose a new improved version where 
the coupling of the f/j-field to dark matter is emphasized at the expense of its self-interaction. This reformulation 
presents the very advantageous possibility that the same potential is used for all galaxies. Using at the same 
time a quasi-isothermal dark matter density and the scalar field helps to better fit the RC of spiral galaxies. In 
addition, new correlations are established. 
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1. Introduction 

Since 1998, many groups have been working on model- 
ing the RC of spiral galaxies with ultra-light scalar fields 
as an alternative to the standard assumption of weakly- 
interacting massive particles. Guzman and Matos (2000) 
presented a model where they assume that the galac- 
tic dark matter consists entirely of a self-interacting real 
scalar field $ endowed with a potential V($) which is 
an exponential function of <f>. First they obtained an ex- 
act static and axially symmetric solution of the Einstein- 
dilaton theory (Guzman & Matos 2000a). This solution 
yields an effective energy density [Idm — V(^) which 
mimics an isothermal dark halo with a constant core den- 
sity. The resulting circular velocity profile of test parti- 
cles involves three free parameters to fit the observations. 
However, the authors (Guzman & Matos 2000b) pointed 
out some difficulties related to the physical origin of the 
kind of exponential potential V($) they have to consider. 
Based on a complex scalar field forming galactic halos af- 
ter Bose condensation, Arbey et al. (Arbey et al. 2001a) 
have suggested a massive (m = 0.4 to 1.6 10 -23 eV/c 2 ) 
and non-interacting complex scalar field, as a dark mat- 
ter candidate alternative to the neutralino, to bind galax- 
ies and flatten the RC of spirals. Their motivation for 
a complex (i. e., charged) scalar field follows from the 
need for stable bounded configurations in the case of self- 
interacting bosons (not interacting with any other kind 



of matter). The model of the authors accounts better for 
the dark matter inside low-luminosity spirals than for the 
brightest objects where baryons dominate. However, be- 
cause they deal with a charged scalar field (Arbey et al. 
2001b), some problems (pointed out by the authors them- 
selves) arise when the limits from galactic dynamics and 
some cosmological constraints are taken simultaneously 
into account. Ureha-Lopez and Liddle (Ureha-Lopez & 
Liddle 2002) have shown that the supermassive black holes 
hosted at the center of galaxies can coexist with scalar field 
halos. Since all the aforementioned studies consider only 
self-interacting scalar fields, the equation satisfied by their 
respective scalar field is of the sourceless Klein- Gordon 
type (including the derivative of the scalar field poten- 
tial). The aim of this work is to show that a real scalar 
field, ip, minimally coupled to gravity sourced by and in- 
teracting with matter (both baryonic and dark matter), 
also may reproduce the RC of spiral galaxies. This can 
be achieved if the force term involved by the scalar field 
can contribute sufficiently to lower the magnitude of the 
orbital momenta of stars and gas, even if ip itself does 
not contribute significantly to the mass of the dark halo. 
The general equation of ip reads (Mbelek & Lachieze-Rey 
2003) 



rv ™ , r dJ . dU 



(1) 



where U = U(tp) and J = J(x a ) denote respectively 
the self-interaction potential and the source term of the 
■0-field ; VyV" is the d'Alembertian of the curved four 
dimensional spacetime endowed with the metric ds 2 = 
g^u dx 11 dx v . The potential U is of a symmetry breaking 
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type but needs not to be specified here ; U{ipoo) = 0, 
f^Ooo) = and §^(-0oc ) > 0, where Voo denotes 
the value of ip at the minimum of the potential U(ip). 
Equation (JTJ describes a minimally coupled classical real 
scalar field with the usual Lagrangian density (Narlikar & 
Padmanabhan 1986) 



£ = 



^ [ i drf> - U - JiP ] y/- det (g^), (2) 



but taking into account that the source term, J, may de- 
pend on the 7/7-field too. This is the case when the coupling 
constants to the other form of matter-energy are general- 
ized to be functions that depend also on the scalar field. 
For our purpose, one may neglect the contribution of the 
electromagnetic fields. Hence, J accounts for the matter 
source only. As usual in the framework of scalar-tensor 
type theories, this is of the general form J = g(ip) T" 
(see e. g., Straumann 2000), where T° is the trace of the 
energy-momentum tensor of the matter source and the 
coupling function g is such that 5(1/' 00) = in the vacuum 
state (minimum of the potential U). 

2. Weak and static field approximation 

Assuming that the excitation Sip = ip — ip^ of the ■0-field 
is small compared to ip^ , the equation of the i/j-field reads 
in the first approximation for a static matter distribution 



ASip - A 



4ttG 

00 5 pi 



(3) 



-1/2 



where A = ^ ^p-(ipoo) J denotes the range of the ip- 
ficld and p is the mass density of the matter fields other 
than the 0-field itself. Out of the matter source, eq© 
reduces to 



ASip - A Sip = 0, 



(4) 



whose solution is the familiar Yukawa potential Sip cx 
exp ( — Xr)/r. Accordingly ip — > ipoo at infinity. Since ip 
is long range (say A > 100 kpc), we may neglect in the 
first approximation the scalar field mass term, A~ 2 , with 
respect to the mass density terms within galaxies. Thus, 
eq0 simplifies as follows 



A / °9 < I W 4nG 



(5) 



1 av 1 > dIlu 1 av — 2 3 ( dx" + dx" dx>- I UB 

fines the Christoffel symbols). In our framework, the 
same procedure applies to an appropriate effective energy- 
momentum tensor T^tj which accounts for the presence of 
the scalar field. Since the trace T£ is involved in the source 
term J, arranging in order to isolate all the contributions 
of the matter in the total Lagrangian density leads to the 
identification T^ f = T» v + g(ip) ip T£ g» v . As a conse- 
quence, this does not lead to the geodesic equation since 
the influence of the scalar field now enters. Indeed, it fol- 
lows for a neutral test body 



ds 



L a/3 a 



d{g{iP)ip) 
ds 



1" - d»{g{iP)iP] 



(6) 



Note that eqJBJ) is generic in scalar-tensor theory like 
the Brans-Dicke (Brans & Dicke 1961) theory or the 
5D compactified Kaluza-Klein theory (Wesson & Ponce 
de Leon 1995). Independently of the physical meaning, 
cqlj^l may be obtained from the variational principle 
= 5 J mc y/g^i/ u v ds with an effective mass m cx 
exp (g(ip) ip). The geodesic equation is recovered in the 
case of a non variable ip-fie\d. Now, in the weak field and 
low velocity limit, one gets Tq ss — \ d^goo (with no 
summation on /1). Hence, in the first order approximation, 
eqlJHJ simplifies to 



dv 
~dt 



= - V [V N 



c 2 ^{iP x )iP Q 



ip 



where the Newtonian potential, Vjv, results from the con- 
tributions of the visible and dark matter components in- 
cluding that of the ip-Held. Newtonian physics is recovered 
in the limit case dip — > 0. 

4. Modelling the RC of spirals 

Since a drag- like force term appears in eqlQ), the rota- 
tional velocity v should read in the general form v = 
\/ Vg + r 2 , where the tangential velocity vg and the radial 
velocity f are derived in spherical coordinates respectively 
from the tangential equation 



1 d(rv g ) dg dip 



and the radial equation 



1 = c m ( r ) 



(r) 



(8) 



(9) 



3. The equation of motion 

The equation of motion of test bodies is obtained following 
the usual procedure (Ciufolini & Wheeler 1995) which con- 
sists of considering a gas of dust (pressureless perfect fluid) 
whose energy- momentum tensor reads T^ v — p^u^u 11 , 
where the density p is a constant. The equation of mo- 
tion is then derived from the conservation of the energy- 
momentum tensor, that is V u T^ v = (where V v denotes 
the covariant derivative V ' V T^ V = d v T^ + T a av T» v + 



where we have set 



™ r) (r) = #f|W>oo)V<oo^ 
ra{r) = m V i S (r) + md a rk( r ) is the total mass up to ra- 
dius r assuming spherical symmetry. In the latter relation, 
Toots (?") denotes the visible matter mass and mdarkir) the 
dark matter mass including both the contributions of the 
dark halo, mhaio{r), and that of the "0-field. So, integrat- 
ing equation (JSJ above yields 



and 



dg 

In (rv g ) = lnj + -^(ipoo) ipoo(ip - tpoo), 



(10) 
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where J would represent the angular momentum per unit 
mass if the ^-field were not excited. Hereafter, we as- 
sume ^(V'oo)V'oo > so that mW')(r) mimics a hidden 
mass profile 1 . In addition, we adopt the model of a quasi- 
isothermal spherical dark halo, that is, a mass density pro- 
file such that 



Phaio{r) = Phaio(ro) ( — ) 2+e 
r 



(11) 



where e is a free parameter which depends on the given 
galaxy and lies within the range ]0, \\. We show below 
that the use of an isothermal dark halo together with 
the possibility of non-circular orbits may explain the ap- 
parent discrepancy with the Navarro- Frenk- White (NFW) 
or Burkert halos assuming circular orbits for the latter. 
Recently, Bournaud et al. (Bournaud et al. 2003) have 
proved for the first time from N— body simulation that a 
large extent of dark matter is necessary for the formation 
of tidal dwarf galaxies. In this study, the authors described 
the dark matter as a standard isothermal sphere. Hence, 
although the main concern of the authors is the tidal in- 
teraction in the outskirts of galaxies, the hypothesis of an 
isothermal (or quasi-isothermal sphere) still proves to be 
worth investigating. For simplicity, we will neglect the con- 
tributions of the bulge and the disk (stellar plus gaseous 
disks) as sources of the ip- field since the dark halo is the 
major mass component of the galaxies. Thus, one finds 
that the static spherical solution ip of eqJSJ reads 



1 dg Gm ha i o (r ) f r 

e oip ro c z 



(12) 



Inserting the solution found for ip in relation l|10|l above, 
the tangential velocity expresses in the form 



ve = vemax G u (x), 



(13) 



where x = r/ro and we have set v = 1/e for convenience 
(see below). The functions G v are defined for x > as 
follows 



G v (x) = — exp — r 

x 



-X/u 



(14) 



These functions are positively defined, bounded from 
above by G„(l) = 1 and satisfy G u (x) k = Gkv{x k ). 
Solving the radial equation ©, on account that r = 
dr/dt = d(^r 2 )/dr (vg = vg(r) and f = f(r) at radius 
r), one finds 



Krf =r(r o y + 2 / ( ^ - G 



m(r) 



(r) 



Hence, one expects m(r) + m^(r) ~ v 2 r/G 
for any good fit to the RC that can be ob- 
tained assuming circular orbits, where 2 (r) = 



(^ffCiAoo) V'ooJ m ha i (r Q ) (r/ro) 1 £ - Actually, the depar- 
ture from circular orbits may help to understand the use 
of the various dark halo density mass profiles, 



Pdm = P* 



U) 



a/3 



1 



0' 



(16) 



that are advocated in the literature to fit the RC of spi- 
ral galaxies, where a = 1 both for NFW and Moore 
09 = i, for the latter) and (3 = 1 both for NFW 
and Burkert (a = 2, for the latter). Indeed, it follows 
from eq® and assuming both conditions r > and 
d(r r 2 /G)/dr > that the quantity rr 2 /G behaves like an 
additional contribution which mimics a mass density pro- 
file as if the orbits were circular (Mbelek 1997). Further, 
m e ff(r) = m(r) + m^(r) + G _1 rr 2 will define an effec- 
tive mass profile consistent with relation i|16|l in the range 
rmin — r ^ r max such that r m i n — ( Phaioir*) / P* ) r* 



and 



(p*/phalo(r*)) r* ; p ha io{r*) < p* 



5. Results 

In the following, we limit our study to those RC for which 
the radial velocity, r, may be neglected with respect to 
the tangential velocity, ve , in the whole range r m i n < r < 
r m ax where we assume that the dark halo mass density 
dominates the baryonic mass density. Furthermore, the 
fits to the individual RC show that v is actually an integer 
n > 2. This is consistent with our former work (Mbelek 
2001) in the sense that combining the solution found for 
ip with relation (|ll|l allows us to express the right-hand 
side of equation @ as the derivative of an effective power- 
law potential which varies as ( if) — f/'oo ) 2 ( n+1 ) . As can be 
seen in figure 1, the greater the n, the steeper is the curve 
y = vg/vemax versus x for x < 1 and the flater it is for 
x > 1. Hence, the steepest is a RC below the turnover 
radius, and it should be flatest beyond. Figure 2 shows 
some fits to individual rotation curves from the samples 
of Rubin et al. (Rubin et al. 1980; Rubin et al. 1985), van 
Albada ct al. (van Albada et al. 1985), Lake et al. (Lake 
et al. 1990). This is achieved by using the least-squares 
fit to search the parameters a and b = In J that yield 
the maximum square, R 2 , of the correlation coefficient for 



dr.(15) 



the relation In (rug) 



a r 



-l/n 



b. Based on the study 



1 Since tp also acts as a stabilizing field (see Mbelek & 
Lachieze-Rey 2003), one may expect that it affects the dy- 
namics within the galaxy by lowering the magnitude of the or- 
bital momentum per unit mass at any given radius. This reads 
rv$<J and consequently &%p < on account of 1101 and by 
assuming ^(t/'oc) ipoo > 0. Consequently, the consistency with 
t/>oo defining the minimum of the potential U(tl>) at r — > oo 
implies that dtp/dr > for Sip decreasing monotonously. 



of a hundred spirals, a is always negative, in agreement 
with Sip < whereas b is always positive. In addition, one 
gets approximately the following statistics : n — 5 for 20% 
of the spirals, 3 < n < 6 for almost a half of them and 
3 < n < 12 for 80% of them. 



2 The case e — 1, which matches a NFW dark halo density 
mass profile far beyond the core radius, yields a point mass 



m 



W0 



(|f (V>oo) ipoo) rn ha i o (r ) that mimics a black hole 



located at the center of the galaxy. 
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Fig. 1. Generic curves y — G n {x) for n — 2,4,8,16,32 
and 64 (from bottom to top). As one can see, the greater 
n, the steeper is the curve below x = 0.5 and the hater it 
is beyond. 





6 8 10 

[arcmin] 



Fig. 2. Rotation curve fits (from top to bottom) for NGC 
4062 (Rubin et al. 1980), NGC 1035 (Rubin et al. 1985), 
NGC 3198 (van Albada et al. 1985; 1 arcmin = 2.68 kpc) 
and DDO 170 (Lake et al. 1990) ; caption : the crosses or 
error bars are the observed data and the solid line is the 
best fit of the RC. 



Figure 3 shows that there is a strong correlation be- 
tween both coefficients a and b. Although displayed here 
for sixteen spirals only, this remains true for the whole 
sample of a hundred spirals we have studied. Expressing 



Fig. 3. Plot of the coefficient b versus the coefficient | a | 
for the galaxies NGC 4062 and UGC 3691 from the sam- 
ple of Rubin et al. (Rubin et al. 1980), NGC 4419, 1035, 
1325, 2742, 3067, 3593, 2639, 4378, 4448, 7606 and UGC 
12810 from the sample of Rubin et al. (Rubin et al. 1985), 
NGC 3198 (van Albada et al. 1985), DDO 170 (Lake et 
al. 1990) and NGC 4051 (Kaneko et al. 1997). One can 
see the strong linear correlation between | a | and b ; the 
crosses are the data from the fits to individual galaxies 
and the solid line is the best fit to the whole sample of a 
hundred galaxies which expresses as b = | a | + 4.25 with 
correlation coeffient R = 0.995. 



r in kpc for all galaxies, one finds 3 b = \ a \ + 4.25 with 
R = 0.995. More generally, there are strong linear corre- 
lations between the integer n and the coefficients a and b. 
Now, the theoretical relations between the fitting param- 
eters (a, b) and the physical ones (ro,V0 m ax) are given by 
r o = (I a I AO" an< l b = n + ln(ro vg max ) . From the lat- 
ter relation, it is clear that the correlation between b and 
n is obtained by the logarithm which alleviates the effect 
of the scatter in the quantity roVg max . As for the corre- 
lation between a and n, it suggests a relation of the kind 

ro=a^l + ^Q , where a ~ 1 kpc and (3 is a dimen- 
sionless positive constant almost independent of the given 
spiral galaxy. Hence, if a natural theoretical explanation 
could be found for the relation between ro and n, only one 
free parameter, namely n, would be needed to fit the RC. 
Since the latter is just an integer, it is tempting to believe 
that this feature may reveal the nature of quantification 
underlying the galactic dark halo mass distributions. 

6. Conclusion 

In this paper, we have shown based on a simple mass 
model that a real scalar field minimally coupled to gravity 
and with matter as its source may reproduce the RC of 
spiral galaxies. This is achieved even if the scalar field 
does not contribute significantly to the galactic dark 



' For the whole sample, Ho is in the range 50 
75 km s" 1 Mpc -1 . 



Mbelek: A model of rotational 



matter mass. Moreover, such a scalar field may also act 
as a stabilizing field. The new correlations that have been 
found open the possibility to fit the RC with only one 
single free parameter which is an integer. Thus, it seems 
possible that long-range scalar fields minimally coupled 
to gravity play a significant role not only at the cosmo- 
logical scale (Mbelek & Lachieze-Rey 2003; del Campo & 
Salgado 2003, Caresia et al. 2004) but also within galaxies. 
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